Symmetry breaking and the random-phase approximation in small quantum dots 



o 
o 

(N 

wo: 
< 



Lloreng Serra, 1 Rashid G. Nazmitdinov, 1, 2 and Antonio Puente 1 

1 Departament de Fisica, Universitat de les Illes Balears, E-07122 Palma de Mallorca, Spain 
liubov Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, 141980 Dubna, Russia 

(Dated: April 10, 2003) 

The random-phase approximation has been used to compute the properties of parabolic two- 
dimensional quantum dots beyond the mean-field approximation. Special emphasis is put on the 
ground state correlation energy, the symmetry restoration and the role of the spurious modes within 
the random-phase approximation. A systematics with the Coulombic interaction strength is pre- 
sented for the 2-electron dot, while for the 6- and 12-electron dots selected cases are discussed. The 
validity of the random-phase approximation is assessed by comparison with available exact results. 
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I. INTRODUCTION 

The random-phase approximation (RPA) is a very suc- 
cessful many-body theory, widely used to describe Fermi 
systems in general. It was originally proposed by Bohm 
and Pines to describe the response function of the elec- 
tron gas and, since then, it has been rederived in many 
different ways and contexts p) . Important developments 
of this approximation for finite systems, including the 
establishment of a formal scheme for both excited and 
ground states, were obtained many years ago in nuclear 
physics Q Unfortunately, great technical difficulties 
have hindered for a long time the solution of the general 
RPA equations, imposing the recourse to severe simpli- 
fying assumptions. In recent times, however, a renewed 
interest in the RPA as a systematic theory of ground 
state correlations in atomic nuclei has appeared (see for 
instance Ref. 0). 

In the electron-gas theory the RPA has played a bench- 
mark role in the development of accurate exchange- 
correlation energy functionals, specially in the high den- 
sity limit 0. When turning to finite electronic systems, 
RPA calculations are generally more involved, due to the 
lack of translational symmetry. Nevertheless, extensive 
applications to metallic clusters (see, e.g., the review by 
Brack and recent studies of molecules have been per- 
formed |6|]. 

A convenient starting point to treat Fermion systems 
is, in many cases, a mean field description like the 
Hartree-Fock (HF) approach. Sclfconsistency between 
the mean field and the single-particle orbitals and total 
energy minimization are the basic conditions at this level. 
It may happen that the selfconsistent solution breaks 
one of the symmetries of the exact Hamiltonian, a well 
known phenomenon called spontaneous symmetry break- 
ing. It has been predicted that the RPA should improve 
the ground state energy. More importantly, it should re- 
store symmetries which could be broken at the mean field 
level. It is our aim in this work to show with a realistic 
calculation for interacting electrons that this is indeed 
the case. We have chosen as an appropriate system a 
two-dimensional parabolic quantum dot, partly because 



of the interest attracted by semiconductor nanostructures 
as a paradigm of artificial quantum systems, and, also, 
because these systems are believed to be good candi- 
dates where electron localization or formation of Wigner 
molecules should be relevant in some regimes (see Ref. 7 
for a recent review). We remind that Wigner molecules 
in quantum dots are associated with the breaking of ro- 
tational symmetry of the mean field. 

The RPA is, in principle, formulated to address the ex- 
cited vibrational states, or vibrons for short, of a Fermion 
system by means of a superposition of particle-hole and 
hole-particle transitions on an a priori unknown corre- 
lated ground state. A basic assumption consists in treat- 
ing each particle-hole pair as a boson, i.e., the so-called 
quasi-boson approximation 0,0,0- Thus defined, the 
RPA spectrum describes the intrinsic vibrations of the 
system and, imposing the condition of excitation vacuum, 
also yields an approximation for the correlated ground 
state. In cases of mean field broken symmetry a remark- 
able property of the RPA is that it separates out the col- 
lective excitation associated with each broken symmetry 
as a zero-energy or spurious RPA mode, with an inertial 
parameter which is also fixed by the approximation. As 
we shall show below, these spurious RPA modes play a 
central role in the restoration of symmetries broken at the 
mean field level. The physical excitations of the system 
consist, therefore, of vibrons and collective modes associ- 
ated with spurious solutions, both kinds being described 
within RPA. 

It should be noted that the restoration of broken sym- 
metries can also be attained via projection techniques 
(see Refs. H,Hj). Examples demonstrating their use for 
the case of quantum dots have been recently presented by 
Yannouleas and Landman |Ioj . In this paper, however, 
we shall focuss our analysis on the RPA description of 
the ground state and its connection to the excited states. 

The electronic islands formed at the GaAs-AlGaAs in- 
terface by using appropriate electrodes constitute a re- 
markable example of controllable quantum system. Many 
properties can be described, in a first approach, by using 
a model of electrons confined by an external potential 
well of parabolic type and restricting the electronic mo- 
tion to the interface plane. These model quantum dots 
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have been intensively studied using a variety of methods, 
including semiclassical, mean field and density-functional 
approaches as well as exact methods attempting a direct 
solution of the Schrodinger equation [?| • Within the lat- 
ter category we find the diagonalization in a basis, the 
coordinate-space and the Monte-Carlo solutions. In the 
exact methods one normally imposes from the beginning 
the Hamiltonian symmetries to the solution and, there- 
fore, the symmetry breaking is not possible (unless de- 
generacy is present). This striking difference with the 
mean-field solution does not necessarily imply that the 
latter one is an artifact. Indeed, the RPA provides a 
consistent scheme for the corresponding physical inter- 
pretation. We attempt in this work a clarification of this 
issue by comparing in a test case, such as small quantum 
dots, exact, mean-field and RPA results. 

The paper is organized as follows: Section II briefly 
introduces the HF mean field approach and discusses the 
appearance of symmetry-broken solutions. Section III is 
devoted to the RPA excitations, introducing the quasi- 
boson approximation and discussing the RPA spurious 
mode connected with a global rotation when the electrons 
are localized in 2-, 6- and 12-electron dots. In Sec. IV we 
discuss the RPA ground state focussing on the correlation 
energy and symmetry restoration, comparing the RPA 
results with exact values. The conclusions are finally 
drawn in Sec. V. 



II. THE MEAN-FIELD SOLUTION 

We consider a system of N electrons whose motion is 
restricted to the xy plane, where a parabolic potential 
induces the formation of an electron island. The full 
Hamiltonian reads 
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where k and m are the dielectric constant and electron's 
effective mass; uq is the external confinement frequency 



and we have used polar coordinates (r 



V 2 )- By 



introducing the length — (li/mwo) 1 ' 2 and energy Eq = 
h 2 /m£Q = Tvusq units we may rewrite Hamiltonian JQl in 
terms of a single adimensional parameter 



R 



as 
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The Rw parameter was introduced in Ref. ^] as a mea- 
sure of the relative importance of electron-electron inter- 
action to confinement potential strength. In Ref. ^3 it 
was used as a measure to define the different phases in 



quantum dots and, in particular, a formation of Wigner 
molecules. In our unit system it is also fulfilled that 
h = m = 1. 



A. The Hartree-Fock equations 

Introducing the basis of oscillator and spin states 
that diagonalizes the first term in Eq. @, { a = 
l,...Af; r\ —1,1} , where a labels the orbital part and 
r\ the spin, an arbitrary single-particle orbital \i) is ex- 
panded as 



(4) 



The oscillator states \a) are characterized by radial (n a ) 
and angular momentum (m a ) quantum numbers -see, 
e.g., Ref.lla We shall assume that each orbital i has non 
zero components only for a given spin orientation rji =f 
or rji =J,, i.e., -Bar) = S^B^. This condition imposes 
good spin symmetry on the single-particle orbitals and it 
is the only symmetry that we shall keep, leaving totally 
unspecified the remaining spatial symmetries. Note that 
the Slater determinant built with these single-particle or- 
bitals will be an eigenstate of the total S z operator but 
not, in general, of S 2 01- ^ n the chosen basis, the HF 
equations are written as a system of nonlinear eigenvalue 
equations for the matrix of B coefficients. The corre- 
sponding equation for the a-th component of orbital i 
reads 



F (0) 



AT 



M / N \ 

£« 

id=i \fc=i / 

,d=l \k=l J . 



. B {i) 
1 ar\i 



(5) 



In these equations is the energy of the oscillator 
state \a) and, analogously, V a bcd and V a bdc are the ma- 
trix elements of the two-body interaction defined above 
(Rw/ r 12) with the corresponding oscillator states. We 
can take advantage of the fact that these matrix elements 
can be obtained analytically in terms of the radial (n a ) 
and angular momentum (m a ) quantum numbers for each 
basis state \a). 

The first and second terms within the square brackets 
of Eq. (|3J are the well-known direct and exchange con- 
tributions to the HF problem, respectively ^5[- We note 
that the direct term is independent of spin, while the 
exchange yields a different contribution for rji =] and 
Vi =[■ Using the selfconsistency approach, both direct 
and exchange terms are determined from the 'preceding' 
iteration and then the eigenvalue problem can be sepa- 
rately solved for each spin subset. Successive repetitions 
of the process lead eventually to the converged solution. 
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It is worth to stress that this formulation of the HF equa- 
tions is particularly well suited to numerical treatment, 
since the use of analytical expressions [l6| for V a bcd min- 
imizes numerical errors. The only remaining source of 
numerical error is, in fact, the precission of the eigen- 
value solver. Nevertheless, a necessary approximation of 
the method is the truncation of the basis size J\f and one 
must check the stability of the results with this parame- 
ter. 



B. Some HF results 

The oscillator basis is defined by including all states 
below a given cut-off E c , i.e., < E c . The value of 
E c has been chosen large enough to satisfy the above 
mentioned convergence criterion. In particular, for N — 
2 and 6 it sufficed to include the first 55 and 153 oscillator 
states, respectively. For N = 12 we used up to 190 states 
which gave good HF convergence. With the purpose to 
present in the following sections a systematic comparison 
of HF, RPA and exact results with varying Rw we have 
taken as test case the N = 2 quantum dot computing its 
properties for Rw = 0.5 to 2.5 in steps of 0.5. Besides, 
some selected cases for N = 6 and 12 have also been 
calculated. 

The HF density 

PHF (r) = ^|^(r)| 2 , (6) 

2, OCC 

where = (r\i) refers to the HF wave functions J3J, 

are shown in Fig. for the N — 2 system. We note 
that between Rw — 1 and 1.5 the system changes from a 
circularly symmetric solution to a spontaneously broken 
one, where the two electrons localize in opposite positions 
in the mean-field frame. These intrinsic localization be- 
comes more and more conspicuous as the Rw parameter 
is increased; an evident manifestation that if electron re- 
pulsion is strong enough the favoured solution consists of 
particles located as far as possible from each other. 

In Fig. |21 we display the density for the N — 6 and 12 
quantum dots, using Rw = 1.89. This value has been 
chosen after Pederiva et al. 17] who, in turn, adjusted 
it to the experiments by Tarucha et ai.JI^. Both cases 
show a clear symmetry breaking; the N — 6 having the 
electrons localized on a ring while the 12-electron dot 
has a central electron dimer surrounded by a ring with 
10 electrons. Regarding the HF total energies, we leave 
the discussion for the following sections where it will be 
compared with RPA and exact results. 

III. RPA EXCITATIONS 

This section presents a brief discussion of the main 
points of the RPA for excitations, which are an essential 
ingredient for the ground state description in the same 
theory. 
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FIG. 1: HF densities for the N = 2 quantum dot with vary- 
ing Rw parameter (shown in the upper left corner of each 
panel). From the outermost contour line inwards each line 
corresponds, respectively, to a density of 0.05, 0.10, 0.15 . . . , 
etc in units of £q 2 . The cutoff in the basis has been chosen 
E c ~ 10.6 Eo, corresponding to a basis size M = 55. 

FIG. 2: Same as Fig. 1 for the N = 6 and 12 quantum dots. 
For clarity the upper plots display a 3d view of the corre- 
sponding densities. The contour lines are defined as in Fig. [5] 
and we have used Rw = 1.89. The basis cutoffs for N = 6 
and 12 are E c = 16 Eq and 18 -Eo, respectively. 

A. The Quasi-Boson approximation 

It is assumed that the system's excitations are created 
by operators of the type 

°l = E a ^ - Y £i a Um) , (7) 

mi 

whose action on the as yet unknown ground state |0) 
yields the excited vibrational states |A) (A > 0); i.e., 

\\)=O{\0). (8) 

Hereafter, we shall refer to the electronic vibrational 
states as vibrons. Using the standard notation we refer 
to particle states (unoccupied HF orbitals) by subscripts 
m, n and to hole states (below the HF Fermi level) by 
i, j. We shall denote the HF Slater determinant for N 
electrons as \v) since it will act as our particle- hole vac- 
uum. The coefficients and in Eq. Q are a 
set of complex amplitudes characterizing each particular 
excitation. 

The quasi-boson approximation (QBA) amounts to 
treat each particle-hole pair (mi) as an elementary boson, 
thus fulfilling 

[a\a m , alaj] w [&mi>&nj] = <>mn% , (9) 

where the first equality introduces the boson operators 
"mi = a ln a i- The result @ is exact in the average, i.e., 
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when considered as expectation value with the uncor- 
rected vacuum \v) as (v\[ala m , a^a^v) — 8 mn 5ij; see 

Refs. 0,0,0. In this approximation, the creation 0\ and 
annihilation 0\ operators obey bosonic commutation re- 
lations 



[Ox,0{,}=S xx , 



(10) 



At the RPA level the eigenstates and eigenvalues are 
obtained by simplifying the exact equations of motion 



H,0{ 



hu x O\ 



(11) 



for the vibron operators O x and excitation energies Tiu>\ 
within the particle-hole space. With the aid of the QBA 

the particle-hole part of the commutator H, 0\ can be 

^Ti ____ l 

derived [2, Li li| and Eq. then transforms to a gener- 
alized eigenvalue problem, well known in the literature. 
We stress that henceforth all commutators are to be un- 
derstood in the above QBA. The RPA equations read 



A B 

B* A* 



= hu>> 



X 

-Y 



(12) 



where the submatrices A and B are given by 



A-mi nj — &mi ^mn^ij 



*mjin *mjni 



-Brni nj — *mnij 



v, 



mnjt 



(13) 



Here, e m i is the energy of each bosonic pair, given in 
terms of the HF eigenvalues as e m i = e m — 

The matrix elements in Eqs. (|13[) may be computed in 
the oscillator basis using expressions similar to those of 
the HF eigenvalue problem. For instance, 



V, 



III J I II 



N 

E T/ r(»™)* rCj)* r(0 r(") 
Vabcd Ba V i B b V] B ->- 

abcd—1 



> cr]i drin 



(14) 



In practice, the numerical solution of Eq. (|12fl requires 
the diagonalization of a rather large matrix, whose di- 
mension depends on the number of hole Nh and particle 
N p states. To be consistent with our imposed spin sym- 
metry we restrict to bosonic pairs where each member of 
the pair has the same spin. 



B. Spurious modes 

It was proved by Thouless 0] that, when the HF 
solution corresponds to a minimum in the energy sur- 
face, the RPA Eq. ltT2)l provides only real frequencies wa- 
Furthcmore, if the full Hamiltonian has any symmetry 
which is broken by the mean field there is a corresponding 
RPA mode with zero frequency (a spurious mode) which 
is orthogonal to the other vibrational (nonzero energy) 



modes. In other words, the generators of symmetries bro- 
ken on the mean field level create eigenstates with zero 
energy in RPA. 

As shown below, we shall be concerned with rotational 
symmetry for which the infinitesimal generator is the 
orbital angular momentum L z . To treat the spurious 
mode related to rotation it is convenient to introduce the 
canonical conjugate operators L z and an angle operator 
$; the latter defined by the following relations 

[H,L Z ] = 0, 
[H,C] = HL Z 
[L Z ,C] = hj . 



operate 

mm 



(15) 

where C — ij§ is an anti-Hcrmitian operator and J is 
the moment of inertia. 

When the above Eqs. are implemented in the RPA 
they lead to the so-called Thouless- Valatin moment of 
inertiajSlJ Jq and to the two RPA vectors 



mi 



■mi "mi ' l mi "mt I j 



a h' ~ r* h 



(16) 



The coefficients of the RPA L z operator are directly given 
by the single-particle HF matrix elements. In contrast, 
for the C operator one needs to solve the linear system 
of equations 



A B \ c 
B* A* { c* 



(17) 



Once these two sets of coefficients are determined, the 
Thouless- Valatin moment of inertia j7o may be calculated 

as 



Jo 



E 

mi 



mi mi 



(18) 



The operators of Eqs. (|l(jfl . together with the RPA 
vibrons 0|, 0\, form a complete set for any operator 
linear in the bosons b mi and b m i. From the point of 
view of its Fermionic character this could represent an 
arbitrary one-particle-one-hole operator of the type 



"■'"mi 



(19) 



The corresponding expansion in terms of RPA excitations 
(for a single spurious mode) reads 



F = E ([Ox,F}0{ + [F,Ol}O x ) 



WA>0 

y Q ([F,c]L z 



[L Z ,F] C) 



(20) 



When applied to the boson operators themselves, the 
above result implies that we can expand each b mi and 
b mi in terms of the vibrons (0\, 0\) and spurious modes 
(L Z ,C). We shall return to this point in more detail in 
the next section. 
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FIG. 3: Evolution of the RPA excitation spectra with the 
number of particle states included. The results correspond to 
N — 2 with Rw = 1 (left) and 2 (right). The same scale has 
been used in both panels for a better comparison. Note that 
increasingly higher RPA excitations appear as N p is increased. 
The same basis size and cutoff of Fig. have been used. 

C. Strength function and sum rules 

The strength function for a general Hermitian single- 
particle operator F reads 

S(E)=J26(E-hu x ) \(MF\0)\ 2 . (21) 

A 

Important quantities related to the strength function are 
its energy moments, usually known as sum-rules. The 
fc-th moment Sk is given by 

S fc =J>WA)* |(A|F|0)| 2 . (22) 

A 

It is well known that the RPA fulfills several exact sum 
rules (see Ref. %2. for recent applications of sum rules to 
quantum dots). In particular, for the dipole operators 
F = Y /l x l or F = Y,i Vi one has Si = h 2 N/(2m). More 
importantly, the RPA also fulfills the Kohn theorem for 
parabolic confinement (see, e.g., Ref. Ufa) stating that the 
only allowed dipole excitation is the center of mass rigid 
motion at an energy fkug. These exact results constitute 
good checks of selfconsistency between HF and RPA, as 
well as of convergence with space size, in numerical cal- 
culations. 



D. Calculated spectra 

Figure displays the eigenvalues hu\ for two cases of 
the N — 2 quantum dot as a function of the number 
of particle states included in the RPA equations. The 
first case (Rw = 1) is circularly symmetric at the mean 
field level while the second one (Rw = 2) has rotational 



broken symmetry (cf. Fig.^). The evolution of the eigen- 
value set with N p shows a remarkable difference in these 
two cases. While all the eigenvalues of the circular dot 
stabilize for a high enough value of N p , the symmetry 
broken solution exhibits one state whose energy keeps 
decreasing as the number of particle states is increased. 
It clearly corresponds to the appearance of an RPA spu- 
rious mode connected with the broken rotational symme- 
try, which in the limit N p — > oo should lie at zero energy. 
Though stable, the first excited state of the circular case 
lies rather low in energy, a hint of the proximity in Rw 
value to the broken symmetry transition point. Actually, 
for Rw = 0.5 the first vibron settles down at a higher en- 
ergy, E « 0.7£ . 

Figure 0] shows the corresponding spectra for the N = 
6 and 12 quantum dots. Both HF solutions correspond to 
broken symmetry cases (cf. Fig. [5| and, therefore, their 
RPA spectra should display spurious solutions. This is 
clearly the situation for N = 6, with a well separated 
low energy mode. For N — 12, however, the situation is 
less convincing and we actually obtain two modes which 
seem to separate from the higher vibron states. We stress 
that for the 12 electron dot the strong increase in num- 
ber of particle-hole pairs prevents us from extending the 
calculation to higher N p values. Nevertheless, a tenta- 
tive interpretation may be given taking into account the 
mean field structure of this dot shown in Fig. [21 One 
would expect a real spurious mode, corresponding to the 
rigid rotation of the foil system, and a soft (low energy) 
mode corresponding to the rotation of the inner electron 
dimer with respect to the outer ring of electrons. These 
two modes could thus explain the behaviour of the two 
lower RPA eigenvalues in the right panel of Fig. 0] 

We analyze next the strength function results. Figure 
|3] displays some representative cases, for the dot with six 
electrons. The upper panels show the electric (x) and 
magnetic dipole (£ z ) strength functions computed within 
the HF method, while the lower ones are the analogous 
results within the RPA. In HF all excitations are above 
the gap between the Fermi level and the first unoccupied 
orbital (the HOMO-LUMO gap) which for N = 6 takes 
the value « 1.8£o- I n RPA, however, one obtains states 
below this gap. Actually, the strength function for the £ z 
operator concentrates almost exclusively on the first RPA 
eigenvalue, at a very small energy; a clear manifestation 
of rotational spurious mode discussed above. For the 



N 


Rw 


Jo (mil) 


2 


1.5 


0.68 


2 


2 


1.26 


2 


2.5 


1.63 


6 


1.89 


16.4 


12 


1.89 


57.6 



TABLE I: Moments of inertia for the spurious rotational 
modes obtained in the Thouless-Valatin approximation (see 
Sec. Iil-B). 
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FIG. 4: Same as Fig. [3]for the N = 6 and 12 quantum dots. 



FIG. 5: HF and RPA excitation cross sections for the N = 6 
quantum dot. Vertical bars show in an arbitrary scale the 
position and height for each state while the solid line displays 
the accumulated contributions to the 5*1 sum rule. Left and 
right panels correspond to the dipole (x) and rotation (£ z ) 
operators, respectively. The inset in the bottom-right panel 
shows an enlarged view, using logarithmic scales, of the RPA 
l z cross section. 



dipole operator there is again a single mode, this time at 
an energy ss 1E , and the total sum rule 51 is ks SE l^. 
As also discussed above, this proves that our results fulfill 
Kohn theorem with an excellent accuracy. 

The Thouless-Valatin moments of inertia obtained for 
the quantum dots with spurious modes are summarized 
in Tab. I. These inertial parameters define physical rota- 
tional bands, characterized by the energies h 2 M 2 / '{2Jq) 
with M = 0, 1, . . . |2^| . We have also checked numerically 
that the Thoules-Valatin moment of inertia coincides, to 
a good accuracy, with the value obtained from a con- 
strained HF calculation for 1Z = H — XL Z as 



Jo 



d 2 {K) _ d(L z ) 



d\ 2 



dX 



(23) 



We stress that the equivalence between the two moments 
of inertia can be fulfilled if and only if a self-consistent 
HF minimum solution is found. 



IV. THE GROUND STATE IN THE RPA 

The RPA ground state |0) is defined from the condition 
that it is the vacuum for all vibrons 



Oa|0) =0 



(24) 



In this section we shall discuss the following aspects of 
the RPA ground state: i) total energy; ii) restoration of 
the HF broken symmetries. 



A. Correlation energy 

The total energy in the RPA (£rpa) can be split in a 
mean field contribution (-Ehf) and a correction (Arpa), 
as 



RPA 



HF 



*RPA 



(25) 



Note that the correction Arpa, with a minus sign, gives 
the standard correlation energy of the system. A de- 
crease in ground state energy (Arpa > 0) will imply an 
improvement with respect to the mean field theory. How- 
ever, since RPA is not based on energy minimization, it is 



FIG. 6: Convergence with number of particle states of the 
RPA correction energy Arpa- The extrapolated line as well 
as the exact value (horizontal line) are also shown. The exact 
values for N = 6 and 12 have been taken from the Monte 
Carlo results of Ref. UtI 



not bound to fulfill the variational principle and, there- 
fore, -Erpa could be even lower than the exact energy. 
In other words, the correlation energy could be overesti- 
mated in RPA, as it has been suggested in the literature 



The RPA correction reads 



A 



RPA 



\ WA>0 / 



(26) 



The above Eq. I|26|l includes the contribution from the 
vibrons at positive frequencies u>\ and, also, from the 
spurious mode associated with a rotational excitation. 
Note that Eq. Il2(j[| is the result of a partial cancellation 
between two large terms. In practice, this may cause 
Arpa to converge rather slowly with space dimension, 
i.e., a high number N p of HF particle orbitals may be 
needed (see below). 



B. Numerical results for ground-state energies 

We have analyzed the RPA ground state for the same 
dots whose excited states have been discussed in the pre- 
ceding sections. Figure displays the evolution of the 
ground state energy -Erpa with number of particle states 
N p . The convergence with N p depends on each specific 
case: while the N = 2, Rw = 1 energy converges fast 
enough, the N = 12 result is clearly still increasing for 
the maximum N p used in the numerics. This is not sur- 
prising, since one naturally expects that the larger Nh 
the larger the space requirements for RPA convergence. 
A similar conclusion may be drawn regarding the inten- 
sity of the interaction, as given by the parameter Ryy- In 
view of the obtained results, -Erpa must be extrapolated, 
in general, from the range of the computationally feasi- 
ble values. We have empirically found that a 4-parameter 
function, a — b/(l + cNp) 1 ^ , leads to quite reasonable fits 
(see Fig. El . 

To assess the quality of -Erpa we compare next with 
available exact energies, i.e., with ab initio solutions of 
the many-body Schrodinger equation. In fact, the stud- 
ied cases were chosen in part to allow for this comparison. 
For the N = 2 dot the exact solution can be obtained 
by introducing center of mass and relative coordinates. 
Only the relative coordinate problem needs to be solved 
numerically and, because of the circular symmetry, it re- 
duces to an straightforward Id radial problem. We have 
solved this equation using a uniform rad ial g rid and stan- 
dard integration methods (see also Refs. l24l . l26l for recent 
discussions of the 2-electron dot exact solution as well as 
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FIG. 7: Results for the 2-electron dot as a function of Rw- 
Right scale corresponds to the occupation numbers of the HF 
hole state (solid symbols) while left scale indicates the ground 
state energies in HF (circles), RPA (squares) and exact solu- 
tion (triangles). See the discussion about occupation numbers 
see in Sec. IV- C 



Ref. [Hj for the application of the projection methods to 
this system). For the larger dots we shall compare with 
the Monte Carlo results of Pederiva et al. [l7| . 

Figure [3 presents a systematics for the N = 2 dots 
with different values of the interaction-confinement ra- 
tio. As anticipated in Sec. IV-A, -Erpa lies below the 
exact value, a manifestation of an overestimated correla- 
tion energy in RPA. Nevertheless, for most of the results 
the RPA clearly improves the HF energy, since -Erpa 
is closer to the exact value than -Ehf- This is not the 
case, however, when Rw = 1; the exact result lying 
almost halfway of RPA and HF in this instance. For 
N = 6 and 12 the results in Fig. increase towards 
the exact values. The extrapolation in these systems 
leads to A RPA (6) = 1.27 E and A RPA (12) = 3.21£ ; 
to be compared with the diffusion Monte Carlo values 
A M c(6) = 1.13.Eo and A M c(12) = 2ME . The RPA 
overestimation of correlation energy is thus of w 12% 
and « 21% for 6 and 12, respectively. The larger differ- 
ence for the 12-electron dot could be partly attributed to 
the extrapolation procedure being based on a relatively 
small space. The deviations with respect to the exact 
values for A^ = 2 arc not surprising since RPA and HF 
approaches are many-body theories aiming at an accu- 
rate description of large enough systems. 
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FIG. 8: Occupation numbers within RPA of the HF single- 
particle orbitals. The step function shown with a dashed line 
separates fully occupied from unoccuppied HF orbitals in the 
mean field picture. 



C. Occupation numbers and expectation value 

In order to construct the RPA ground state we use Eq. 
(|24H and the condition 

£*|0)=0, (27) 

that insures the rotational invariance of the ground state. 
One seeks solutions of the form 0, H, H 

|0) = N a e 5 \v) (28) 

with the S operator involving the creation of two bosons, 

S =l E Zminjbl^. (29) 
mi nj 

In Eqs. (|28l) and (|29|l N is a normalization constant and 
the matrix Z min j is, in general, complex and symmetric 
in the boson indexes, i.e., Z m i n j = Z n j m i. 
Using the general indentity 

Fe s = e s (f+[F,S] + ±[[F,S],S] + ..'\ (30) 



for the operators F = 0\ and F = L z ; as well as Eq. 
Eq. 1(23) and the QBA one finds 

(F+[F t S])\v) = 0. (31) 

For the case of non-broken symmetry the above require- 
ment leads to the following equations for the Z m i n j co- 
efficients 

V (X)* _ y y-(A)* /on\ 

nj 

When the mean-field breaks rotational symmetry we have 
to complement the above system of equations for vibrons, 
which are reduced by one equation, with the additional 
condition for the spurious mode 

nj 

Our interest is focused on the occupation numbers 
of the HF orbitals in the correlated ground state, 
(Q|c4a m |0) and <0|aja,|0). Using 

Wl^bU = Sijbl (34) 
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one obtains 



(014^10) - E^-ni(o|6L^|o> 

inj 

<0|oto,|0) = 1 - J2 Z minj (OI&L^IO) • (35) 

mnj 

By substituting in the above expressions the expansion of 
the boson operators 6^ in terms of vibrons and spurious 
modes, Eq. I|2U|). we finally obtain 

(0|aLa m |0) = ^£|^ + (0|$ 2 |0>^K«| 2 ( 36 ) 

iX i 

<0|a.k|0) = 



^A |2 
1 mi I 



(0|$ 2 |0}^K, 



, ,2 

/// / I 



(37) 



These equations generalize the result known in the lit- 
erature [2j by introducing an additional term related to 
the canonical variables of the spurious mode {L z , $}. It 
should be noted that the factor 1/2 in Eqs. (|36|l and 
(I37|) is introduced following Ref. |22l where the occupa- 
tion numbers were calculated using Fermionic anticom- 
mutator rules, without referring to the QBA. 

The above discussion can be easily extended to obtain 
the expectation value of any 1-body operator such as, 
e.g., the particle density 



a/3 



(38) 



where indexes a and (3 run over all the HF set of orbitals 
and p Q ^(r) = tp* a (r)( y 9 / g(r) (with Lp a the HF wave func- 
tions). Omitting for the sake of presentation the spatial 
dependence, one can write the RPA ground state density 
as 



(§\p\$) = ^ Pli - Ap h + Ap v 



(39) 



|20j. Therefore, we propose the following procedure to 
fix its value. Treating the unknown matrix element as a 
parameter (7 = (0|$ 2 |0)) we introduce the function 



F(i) = J dr [p( 7 ,r)-p( 7 ,r)r (42) 
where ,0(7, r) refers to the RPA density Eq. (|39f) and 

/5(7,r) = i/ <0p(7,r) (43) 
27r Jo 

is the corrresponding angular average. Since |0) is ro- 
tationally invariant, it is natural to require the physical 
value 7 = 70 to be a minimum of F(-y). Therefore, it 
fulfills the condition 



= . 



(44) 



7=70 



Resolving this equation we obtain the unknown quantity 
as 



7o 



A 
B 



(45) 



where 
A = 

B = 



r 22 

/ dr^^(-l) fe+m K-a fc )(& m -& m ) 

« 7 n ™ — 1 



fe=0m=l 
2 



f drl Y,(-l) m (b m -b rn ) 

\m=l j 



(46) 



and we have defined the a and b coefficients according to 
Eqs. (gOJ and (HU as 



Ap h = m + 761 
Ap p ee a 2 + 7&a 



(47) 



with 



Ap h = 



Ap p 



(A)* 
mi 



Y, 



mj 



Xijm 



<0|$ 2 |0> 



ijm 



'mi 



V Y W Y {x) * 

Xmnj 



- mj nj 



(0|$ 2 |0) 



^ ^ Pmn 



{z)* p(z) 
mj c nj 



(40) 



(41) 



The occupation numbers Eqs. I|36|l and l|37(l . and den- 
sity expectation value Eq. (|39|l depend on the matrix 
element (0|$ 2 |0). Involving the two-body operator $ 2 
this matrix element is not determined within the QBA 



The tilde coefficients are defined through analogous equa- 
tions with the circularly averaged densities. 



D. Numerical expectation values 

The above discussed Fig. [7| contains also the numeri- 
cal values of the hole state occupation for N — 2 -there 
is only one occuppied orbital (spin degenerate) in this 
case. The occupation is close to 1 for weak interaction 
(low Rw) and it decreases in an important way as the 
interaction strength is increased. This proofs that the 
RPA ground state deviates from the HF one when in- 
teraction becomes strong and correlations allow particles 
and holes to be found above and below the Fermi level, 
respectively. The occupations for the other two dots an- 
alyzed in this work are shown in Fig. |SJ In these cases 
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FIG. 9: Comparison of RPA and exact radial densities for the 
N — 2 dot with different Rw parameters. In the Rw = 0.5 
and 1 cases the HF density, which is circular, is also displayed. 



we have two (for N = 6) and four (for N — 12) different 
hole levels with occupations varying with their respective 
energies. In general, the RPA occupation numbers show 
a smoother distribution than the HF ones, the maximum 
variations occurring nearby the Fermi level. 

The RPA density is compared with the exact one in 
Fig. for the 2-electron dot. For all Rw parameter val- 
ues the RPA density is circularly symmetric; fulfilling the 
symmetry restoration discussed above. Therefore, in Fig. 
we focus exclusively on the radial dependence. While 
the dot edge is well reproduced, a conspicuous feature is 
the underestimation of the central density by the RPA, 
specially at large Rw values. Our result is in qualita- 
tive agreement with that of Reinhard [2^ for electrons 
in jellium spheres (where a circularly symmetric mean 
field was imposed) . As shown by Reinhard the RPA cor- 
rectly describes the low q components (large r's) of the 
form factor F(q) — J dr exp (iqr)p(r), but it fails for the 
large q contributions (low r's). The overcorrection of 
both central density and correlation energy seem thus to 
be peculiarities of the RPA. 

Analogous results for the 6- and 1 2-electron dots are 
presented in Fig. EH An excellent restoration of the cir- 
cular symmetry is obtained for the N — 6 dot -we remind 
the reader that the corresponding HF density was shown 
in Fig. For N = 12 the RPA density, though more 
circular than the HF one, still has some residual defor- 
mation. This can be surely attributed to incompleteness 
of the RPA space considered in the numerical calculation. 
In fact, for this dot the spurious mode separation (Fig.|3J, 
the ground state energy (Fig. [SJ , and the density (Fig. 
llOf) are all indicating that convergence with number of 
particle states is rather slow and, therefore, very difficult 
to be achieved in a numerical calculation. We emphasize 
that circular symmetry can be restored within the RPA 
only if the contribution from the spurious mode is taken 
into account. 



FIG. 10: RPA symmetry restoration of the HF densities dis- 
played in Fig. |21 



V. SUMMARY AND OUTLOOK 

The RPA is a general framework where mean field 
theory can be improved. Focussing on two-dimensional 
quantum dots we have reviewed the description of excited 
and ground states within the RPA. Some of the results 
known in the literature have been generalized to include 
broken symmetries at the mean field level. Quite im- 
portantly, the RPA provides a physical interpretation of 
this broken-symmetry phenomenon: mean-field solutions 
lacking one symmetry of the Hamiltonian represent, ac- 
cording to the RPA, true internal deformations of the sys- 
tem having an associated collective motion at vanishing 
energy (the spurious RPA mode). The ground state en- 
ergy, occupation numbers, as well as general expectation 
values with 1-body operators, including the spurious- 
mode contribution, have been obtained. We proposed 
the procedure to restore the rotational symmetry bro- 
ken at the mean field, which can be extended for other 
symmetry breaking cases. 

Small quantum dots have provided us a good scenario 
to assess the above properties of the RPA by numerical 
calculations. For large enough values of the interaction- 
confinement ratio Rw the HF mean field breaks circular 
symmetry; the electrons being localized in specific geo- 
metric distributions. In these cases we obtain an spurious 
RPA mode related to the global rotation, with its asso- 
ciated moment of inertia. The corresponding rotational 
spectra may be associated with the rotation of a Wigner 
molecule |2||. A systematics with Rw for the 2-electron 
dot has been presented, while for 6- and 12-electron dots 
the calculations have been performed with an Rw value 
suggested by the experiments. 

When comparing with exact results, the RPA gener- 
ally overestimates the correlation energy and the central 
dot densities by an amount which depends on the elec- 
tron number and the Hamiltonian parameters. Another 
general behaviour of RPA is the slow convergence with 
space dimensions, requiring in practical applications a 
very high number of empty HF orbitals. For systems 
with many electrons this renders the numerical calcula- 
tion quite a formidable task. In spite of these difficulties 
the RPA corrections essentially improve the HF energies 
with respect to exact values. It should be noted that 
the QBA only approximately satisfies the Pauli principle 
(see Eq.@). To solve this deficiency of the QBA could be 
one way to improve the description of the ground state. 
On the other hand, a more systematic treatment of the 
ground state energies would require, probably, alterna- 
tive many body techniques. 

The exact spin symmetry has been imposed on the 
single-particle states in this work, both at the HF and 
RPA levels. An interesting extension is, therefore, to re- 
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lax this constraint allowing the HF orbitals to be general 
two-component spinors. The allowed symmetry break- 
ing would eventually lead to ground states consisting of 
vectorial spin textures, lacking a single quantization axis. 
This type of states have been predicted for quantum dots 
in large magnetic fields and the question arises what col- 
lective spurious modes will be associated with them. This 
and related issues are left for a future work. 
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